Gaussian superconducting fluctuations, thermal transport, and the Nernst effect 
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We calculate the contribution of superconducting fluctuations to thermal transport in the normal 
state, at low magnetic fields. We do so in the Gaussian approximation to their critical dynamics 
which is also the Aslamazov-Larkin approximation in the microscopies. Our results for the ther- 
mal conductivity tensor and the transverse thermoelectric response are new. The latter compare 
favorably with the data of Ong and collaborators on the Nernst effect in the cuprates. 



The study of fluctuations in superconductors 0] had a 
marked revival after the discovery of the cuprate super- 
conductors, ft was realized early that their short coher- 
ence lengths produce a large regime of strong fluctuation 
phenomena Q] • More recently the "pseudogap" region of 
the cuprate phase diagram has come into sharp focus, and 
one line of thought attributes its features to strong super- 
conducting fluctuations 0, ■ Among the most striking 
experimental findings that plausibly support this link are 
recent measurements by Ong's group [5( where the Nernst 
effect becomes sizeable far above T c in the underdoped 
regime of the high-temperature superconductors. As the 
Nernst effect is small in ordinary metals and large in the 
vortex state of superconductors the a priori case for cred- 
iting superconducting fluctuations is strong. 

In this Letter we calculate the contribution of Gaus- 
sian superconducting fluctuations to the transverse ther- 
moelectric response above T c in the low magnetic field 
limit — the simplest computation that can be used for 
a quantitative comparison with these experimental re- 
sults. This is already a non-trivial exercise as it requires 
a proper subtraction of bulk magnetization currents that 
enter naive expressions [(J. We also give results for the 
thermal conductivity tensor, thus presenting a full pic- 
ture for Gaussian thermal transport in the system. 

Our result for the transverse thermoelectric response 
depends only on the superconducting coherence length £, 
making it particularly suitable for comparison with ex- 
periment. We present a comparison with data from dif- 
ferent samples of La2_ 2; Sr 2 ;Cu04 (LSCO). We find that 
the Gaussian Nernst effect, using the actual T c , is quan- 
titatively comparable to the measured signal in the opti- 
mally doped and overdoped samples. For the underdoped 
sample the measured signal in the "pseudogap" region is 
larger and requires postulating a suppression of the ac- 
tual T c from the mean-field T c , consistent with the su- 
perconducting fluctuations scenario, in order to achieve 
an understanding of its magnitude, as we illustrate by a 
self-consistent Hartree computation. 

Before proceeding, we note that the traditional de- 
scription of the Nernst effect is of thermally driven vor- 
tices producing a transverse voltage via phase slips. In 
a fluctuation regime this description is only convenient if 
the vortices provide an effective parameterization, as we 



expect will be the case near a two dimensional Kosterlitz- 
Thouless transition. At higher temperatures, where the 
vortices and anti- vortices are strongly overlapping, a dif- 
ferent simplification involving Gaussian fluctuations be- 
comes available and we will take recourse to that in the 
present paper. 

Gaussian computations. — Our calculations can be car- 
ried out in two equivalent ways. First, they involve fol- 
lowing Aslamazov and Larkin [?J and keeping the Feyn- 
man diagrams that now bear their names and arise in a 
Gaussian treatment of the quantum functional integral. 
Second, as was noted, e.g., in the case of paraconductiv- 
ity 0, one can keep Gaussian fluctuations in a stochas- 
tic time-dependent Ginzburg-Landau equation (TDGL) 
that is designed to recover the equilibrium Ginzburg- 
Landau free energy. At this level, with the coefficients 
of the TDGL derived from BCS microscopies near the 
mean-field T c , the two computations are identical. The 
TDGL, however, has a second interpretation — as a model 
of the critical dynamics traditionally assumed for super- 
conductors (model A) || which should have a wider va- 
lidity, with the coefficients no longer constrained by BCS 
microscopies. In the following we will mostly use the 
TDGL description, while appealing to the microscopies 
to justify the form of the current operators |9j. 

The stochastic TDGL is defined by the Ginzburg- 
Landau free energy 
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where a — ao(T — T c ), and the time evolution 
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Equation J2J describes the decay of the order parame- 
ter configuration towards its minimal free energy, inter- 
rupted by thermal fluctuations introduced through the 
white noise C with correlator 



(C(r,t)ar',t'))=2Tr6(r-r')8(t-t'). 



(3) 



The relaxation time of the order parameter is r, and r' 
leads to the precession of the order parameter during the 
relaxation process (in this convention, t/H and t'/H are 



2 



dimensionless) . In the presence of particle-hole symme- 
try t' = 0. 

To calculate the response of the system, Eq. is 
to be solved in the presence of the appropriate driving 
field (electric field or temperature gradient). Note that 
a temperature gradient appears in the TDGL both in 
the parameter a and in the noise correlator, Eq. At 
the level of Gaussian fluctuations the TDGL is a lin- 
ear equation, with a solution of the the form tp(r,t) — 
J dr' db' G{r,t;r' ,t')((r' ,t'), where G is the Green func- 
tion of the TDGL. This expression may then be used to 
calculate the electric and heat currents, 
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Ullah and Dorsey [lfj used precisely the same model. 

Currents and fields. — As model A has no conservation 
laws, the traditional route to identifying currents is not 
available. While it is possible to construct hydrodvnarnic 
arguments that justify the forms in (@J and (0 1 1| , for 
our purposes it is sufficient to appeal to the microscopies 
of the Aslamazov-Larkin contribution, and show how the 
heat current arises from the appropriate vertex in the 
microscopic theory. [l2T | 

The microscopic theory is perhaps best recast in a 
functional integral approach. By means of the Hubbard- 
Stratonovich decoupling of the interaction, the expecta- 
tion value of a current operator is expressed as an imag- 
inary time functional integral over the pairing field, 



0) 



-S off (i/>») 



/ Dtp Dtp e-SoffW-.V^) 



(6) 



where S e s(ip, ip, </>) is the effective action for the pairing 
field t/}, in presence of a potential <f>, and ( j )j 1 $d> is the 
current of a non-interacting Fermi gas driven by elec- 
tric and pairing fields. Ignoring the dependence of ( j ) 
is the Aslamazov-Larkin approximation; the expectation 
value is the vertex in their diagram and its long wave- 
length, low frequency form is the TDGL current. For 
the electric current this was done in the original calcu- 
lation 0- Using the microscopic heat current operator, 
one obtains a relation between heat current and electric 
current vertices, J Q = — wJ e /2e, independent of disorder. 
This corrects a factor of 2 in the calculation of Reizer and 
Sergeev E! The heat current in the TDGL, Eq. ©, is 
now immediate by this result. 

To obtain a response to a temperature gradient in 
the microscopic theory, one calculates the response to a 
"gravitational field" , coupled to the energy density in the 
Hamiltonian, and then uses an Einstein relation to ob- 
tain the response to a temperature gradient |l4| . The mi- 
croscopic calculation of the Aslamazov-Larkin diagrams 



TABLE I: Contribution of Gaussian fluctuations to ther- 
moelectric and thermal transport coefficients. Here, l B = 
(hc/eB) 1 / 2 is the magnetic length, f = h/(2m* a) 1//2 is the 
Ginzburg-Landau coherence length, T A is a cutoff tempera- 
ture, and e* = — 2e is used. Exact coefficients are given for 
diverging terms only (for non-diverging contributions, singu- 
lar behavior is indicated without specifying prefactors). 



two dimensions 



three dimensions 



sc J_ £ ll f t e f 

axy fah7[ cx t~¥ c 12^ ft 4 K y/r=Tc 



k s x ° co -ci(T-T c )ln 



T-T £ 
f eB hr' , / Tj 



- a(T-T c ) 3 / 2 



K xy ~. ; 5" In 

47r m*c t 1 



"Previously considered in Refs. |l(lll.'llTfjl . 

leads to the same results obtained using the TDGL ap- 
proach, Eqs. 0-©, where the response to a tempera- 
ture gradient is obtained directly. Ignored here are the 
Maki-Thompson and den sity of state corrections to the 
normal state response [l|, lf5j . 

The transport coefficients are defined by the standard 
linear response relations, 
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The off-diagonal thermoelectric tensors obey the Onsager 
relations a = Ta. Below, we calculate the contribution 
of Gaussian superconducting fluctuations to the thermo- 
electric and the thermal conductivity tensors (a sc and 
k sc ). We do not consider here the electrical conductivity 
tensor er sc , which has been studied extensively 0. 

Before proceeding, we note that particle-hole symme- 
try (t' = 0) implies that = a s x c x = n X y = 0. In 
the following wc break particle- hole symmetry only when 
considering these coefficients. In addition, wc calculate 
the longitudinal coefficients in the absence of magnetic 
field B, and the transverse coefficients to linear order in 
B (valid for B^ 2 <C hc/e*). The final results in both two 
and three dimensions are presented in Table [I] The gen- 
eralization of these results to a layered superconductor 
using a Lawrence-Doniach model is straightforward. 

The thermoelectric response a% c x may be calculated in 
two ways, either as the heat current response to an elec- 
tric field or as the electric current response to a tempera- 
ture gradient. In this way Onsager relations are verified. 
The result has a logarithmic divergence at T c in two di- 
mensions on ei Ei 

The calculation of the transverse thermoelectric coef- 
ficient a x ^ raises the issue of equilibrium magnetization 
currents. To illustrate this point, consider the result of 
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calculating both the heat current response to an electric 
field and the electric current response to a temperature 
gradient (in two dimensions), 
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Equation agrees with the calculations of Ullah and 
Dorsey |irj|. However, the two results, Eqs. 0-©, give 
a different answer for the transverse thermoelectric re- 
sponse, and are clearly incompatible with Onsager rela- 
tions. The reason for this apparent discrepancy is that 
Eqs. give the total currents in the system, con- 

taining contributions from both transport and magneti- 
zation currents. In the context of superconductivity this 

issue was raised in Refs. nana, 

but was later ignored 
in works regarding fluctuation contributions Jiol fig . |20| . 
Here we follow Cooper, Halperin, and Ruzin 6] who con- 
sidered the general problem of magnetization currents in 
magneto-thermoelectric transport in detail (they encoun- 
tered it in the context of the quantum Hall effect). 

In the presence of a magnetic field, the system has 
magnetization currents in equilibrium. These currents 
are divergence-free, and as a consequence do not make 
any contribution to the net current flows that are mea- 
sured in a transport experiment. The total currents [cal- 
culated by Eqs. 0-0] are thus a sum of transport and 
magnetization parts, 
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In the present case , following the arguments of Ref. p , 
the magnetization electric and heat currents which con- 
tribute to the total currents in Eqs. 0-0 are 
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where M is the equilibrium magnetization. The fluctua- 
tion contribution to the magnetization is found either by 
thermodynamics or by calculating the current flowing 
in the system in equilibrium in the presence of a varying 
system parameter (such as the magnetic field or T c ). 

By subtracting the magnetization currents, Eq. 
we obtain the result for a%° (see Table QJ, which diverges 
as the conductivity, a%° oc cr§° oc (T — T c )( d " 4 )/ 2 . As ex- 
pected, Onsager relations are recovered by this subtrac- 
tion. We note that the correction due to magnetization 
currents is not small: In response to an electric field, the 
magnetization current accounts for two thirds of the total 
heat current in the weak magnetic field limit. In response 
to a temperature gradient, the magnetization current is 
more singular than the transport current at T c . 

Finally, we consider the thermal conductivity tensor 
k. The longitudinal thermal conductivity k£° does not 
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FIG. 1: Points are a xx (u — u n ) for different samples of 
La 2 - ic Sr ir Cu04 [13, with x = 0.12 (underdoped, T c = 29 K), 
x — 0.17 (near optimal doping, T c = 36 K), and x = 0.2 (over- 
doped, T c = 27 K). Solid line is theoretical value of a X y/B, 

using = 30 A and an anisotropy of 7 = 20. Dashed line 
is obtained using a Hartree approximation (see text). 

diverge in any dimension, although it is singular. (In 
Ref. |2l| a diverging result was obtained due to the wrong 
form of the heat current vertex.) The transverse thermal 
conductivity has a logarithmic divergence at T c in 
two dimensions (similar to a s x c x ). In contrast with an 
electric field, a temperature gradient does not give rise 
to a magnetization heat current in the present case. 

Cuprate data. — We now compare our results for a x ^ 
with experiment. Consider the measurement of the 
Nernst effect: The sample is placed in a temperature 
gradient (— VT) || x in the presence of a magnetic field 
B || z. The transverse electric field E y is measured in the 
absence of any transport electric current. Imposing the 
condition j£ r = in Eq. gives Nernst signal, 
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Taking the values for a and a from the Gaussian contri- 
butions, we find that in this approximation the Nernst 
effect tends to a constant at T c , v(T c ) = a|°/cr|°B. How- 
ever, for comparison with experiment the fluctuation con- 
tributions, cr sc and a sc , should be added to the normal 
state contribution, a n and a n . Moreover, we will con- 
sider temperatures which are not too close to T Cl such 
that a sc <C cr" w a. Equation l|12|) may then be written 
in the low magnetic field limit as 
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(13) 



The normal state Nernst effect v n is generally small due 
to a cancellation between the two terms in Eq. i|12|) . 

Eq. I|13fl is now used for a comparison with experiment 
in the low magnetic field limit, presented in Fig. ^ f° r 
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data measured on three samples of LSCO [22j . The ex- 
perimental data shows a xx {v — v n ), where v and a xx was 
measured at each temperature, while v n is deduced from 
an extrapolation from the high temperature regime. As 
v n is very small in this material, this extrapolation is 
expected to introduce a relatively small error. 

On the theoretical side, the result for a x c y (see Ta- 
ble HJ) is generalized for a layered superconductor using 
the Lawrence-Doniach model, 
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Here, s is the interlayer spacing, and £ Q f, (£ c ) is the co- 
herence length in the direction parallel (perpendicular) 
to the planes. Note that the only fitting parameters are 



the coherence length ^ (£ ab = ^>y/T c /(T - T c ) ) and 
the anisotropy 7 = £ a b/£c- In Fig. {Qa^/B is plotted 

using ^ = 30 A and 7 = 20 (cf., e.g., Ref. 0). 

As this comparison suggests, Gaussian superconduct- 
ing fluctuations are sufficient to explain the magnitude 
of the observed signal above T c in the low magnetic field 
limit in the optimally doped and overdoped samples. 
Two questions arise immediately: what is special about 
the cuprates and why is the fluctuation contribution to 
the Nernst signal dominant at such high temperatures? 
From our results, the answer to the the first appears to 
lie in their anisotropy and in their smaller conductivity 
near T c which together boost v by two orders of mag- 
nitude from values we would predict for bulk Al or Nb, 
despite their larger coherence lengths. In this context 
it would be interesting to study low temperature super- 
conductors with high resistivities, such as the Nb films 
studied for fluctuation effects by Hsu and Kapitulnik [24| , 
where our formula would indicate a Nernst signal compa- 
rable to optimally doped LSCO. In answer to the second 
question we note that the Nernst signal is particularly 
suitable for a fluctuation measurement due to the small 
background signal in the normal state. 

In the underdoped sample the measured signal is larger 
and cannot be explained by Gaussian fluctuations if one 
uses the actual T c . However, this is the pseudogap region 
and in the interpretation considered in this paper, the 
actual T c is suppressed from the mean-field transition 
temperature T™ F by non-Gaussian fluctuations and so a 
naive fit is not justified; indeed, there is presumably a 
small effect of this kind even at optimal doping which is 
unimportant in the region considered in our fit. 

The increase in fluctuations with underdoping can 
be modelled by a growing quartic term in the 
Ginzburg-Landau functional as well as by increased two- 
dimensionality while the increase in T C MF needs to be put 
in directly. To get a sense of what these would do, we 
treat the resulting problem via a self-consistent Hartree 
approximation (see, e.g., Ref. 0), which should be valid 
to some extent below T" F . This amounts to replacing a 



by the self-consistent solution of a = a + b(\ip\ 2 ). That 
this is plausible is demonstrated in Fig.^ where we fit the 
data for the underdoped sample (using a two-dimensional 
Hartree approximation, with T C MF = 49 K, T A = 200 K, 
b/a = 0.63h 2 /m*, and = 30 A). 
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